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Abstract 



We use three-point function QCD sum rules to calculate the form factor 
Aq(o 2 ) appearing in the matrix element of the flavour-changing axial current 
between the D{B) state and a vector meson state. We describe the role of this 
form factor in nonleptonic D(B) — > PV decays and analyze the light SU(3)f 
symmetry breaking effects. We also discuss a proposal to relate the branching 
ratio of B — > K*^ to the spectrum of the semileptonic B — ► plv decay. 
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I. INTRODUCTION 



The theoretical description of weak exclusive nonleptonic decays of heavy mesons is 
carried out, in the factorization approximation, in three steps. First, an effective hamiltonian 
is constructed taking into account the effects of hard gluon exchanges |JJ. Then, the hadronic 
matrix elements are factorized into the product of current-particle matrix elements that can 
be either inferred from experiment or calculated theoretically (notice that we only consider 
two body decays) |2|||. Finally, strong rescattering effects (mainly for D meson decays) 
are included considering the coupling to intermediate resonances |3j] or using the measured 
phase shifts 0. 

The second step is dictated by our inability in reliably computing matrix elements of 
the effective hamiltonian between the external states involved in the decay. Let us con- 
sider, for example, the transition D + (p) — > K (p')^ + (q), which is governed by the effective 
hamiltonian: 
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= 7^ V ud V cs 



^o 1 + ^^o 2 



where Gp is the Fermi constant, Vhk are CKM matrix elements and the local operators 0\ 
and O2 are given by: 

Or = s 7m (1 - 7s )c wf(l ~ J 5 )d (1.2) 
2 = S7„(l - l5 )d u 7 "(l - 7s )c • (1.3) 

The Wilson coefficients c + (/i) and c_(/x) account for the hard gluon effects in the renormal- 
ization of Hw from my/ to the low energy scale \i\ they are known at the next-to-leading 
order approximation ||. The amplitude A(D + — ► K ix + ) is written, in the factorization 
approximation, as follows: 

A{D + K*\+) fact = -9£v* d V cs < T\ + \ [c.O, + c 2 2 ] \D + > fact 

= -ffOM (ci + §-)< K*°\s lfl (l - 75 )c| J D+ >< n+\uY(l ~ ls)d\0 > 

+ (ca + ^) < ir + \u^(l - 75 )c| J D+ X K*°\s^(l - 75 )d|0 > } (1.4) 
where c\ = (c+ + c_)/2, c 2 = (c+ — c_)/2. To obtain ( |1.4|) we have used the identities 
Ox = ^0 2 + l - [»7 M (1 - l5 )X a d] [uY(l - 75 )A a c] 

O2 = + \ [*y„(l - 75 )A a c] [uY(l - l,)X a d] (1.5) 

where N c is the number of colours and A a are the Gell-Mann SU (3) c matrices. Therefore, 
the amplitude is expressed in terms of current-particle matrix elements 
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< K*\p', A)|s 7M d|0 >= f K *m K ^{\) (1.6) 

<n + (q)\u-f^ 5 d\0>= -if n q^ (1.7) 

(e(A) is the if* polarization vector) and of the matrix elements governing the semileptonic 
transition D — > K*{i:)tv: 



< T°(p', A)^ 7m (1 - 75 )c| J D + (p) >= {w^V^^^ 

^ 777. 7~t — k T 



2_ 



> D + m x 0^i(9 2 )< - ^ g2) (e* • P )(p + p% 

m D + m K * 



- (e* • ^^g^g 2 ) - A (q 2 ))}} , (1.8) 

< 7r + (q)\uj,c\D + (p) >= (p + q)^' 2 ) + ^^[F (p' 2 ) - ^(p' 2 )^ (1.9) 

(q = p — p'). Notice that we have used the BSW |7j parameterization of the form factors; in 
this parameterization A\, A 2 and A 3 are not independent, but they fulfil the relation: 

2 m D + m K , 2 m D - m K * 2 

^3? = — = = A 2 {q 1.10 

2m k* 2mx* 

with 

A 3 (0) = A (0) (1.11) 

in order to avoid the unphysical singularity at q 2 = in eq. (|1.8|); for the same reason 
F (0) = F 1 (0). 

The amplitude of the process we are considering is therefore reduced to: 
A(D+^K*\+) fact = -^V: d V cs 2rn K *(e*.p) { ( Cl + ^)/ 7r A (m 2 ) 

Cl 



+ (c 2 + ^)/K^i(m|.)}- (1.12) 



The leptonic constants and fx* can be derived from experiment. As for the Wilson 
coefficients, it has been observed fHH that consistently neglecting all 1/N C terms in the 
factorized amplitudes like eq. ( pL . 1 2| ) , and using the form factors e.g. from the BSW model 



0, the overall agreement with the available experimental data seems to improve. This is 
an appealing feature of this approach, since factorization becomes exact in the multicolor 
chromodynamics, in the limit N c —>■ oo ||. Dynamical justifications of the rule of discarding 
1/N C contributions have been investigated by analyzing the sign and the size of nonfactor- 
izable matrix elements in some D and B meson nonleptonic decays [H]|ll[ . 
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However, in order to have a quantitative overview of the validity of the factorization 
approach we need a careful determination, either from experiment or by a QCD calculation, 
of the form factors appearing in the factorized amplitudes like ( [L.12| ). This is an important 
preliminary study in the analysis of the decay channels where naive factorization seems to 
fail |2] . 

Early estimates of F\, V, A\ and A 2 , evaluated at q 2 = or at q 2 = q^ax, are avail- 
able from constituent quark models |7| JTB|JT^| . More recently, these form factors have been 
computed, at a fixed value of q 2 , by QCD sum rules [|15H T7jl and by lattice QCD fl8|| ; they 
have also been derived in the framework of models that incorporate heavy quark and chiral 
symmetries Moreover, QCD sum rules have been used to determine the q 2 dependence 
of Fi, V, Ax, A 2 and F |l5Hr^, pO| , p]] ; data on this functional dependence are also available 



from lattice QCD, but with large error bars [TLB 



As for A (q 2 ), the procedure adopted so far consists in using eqs. (|1.10[ - |l.ll|) to get A (0) 
from ^4i(0) and ^(O), and then in assuming a suitable dependence on q 2 , invoking the 
dominance of the nearest pole in the t = q 2 channel. 

This procedure has two difficulties. The first one is that Ax(0) and A 2 (0) are predicted 
by QCD sum rules and lattice QCD in a range of values; the uncertainty (larger for A 2 ) is 
determined by the variation of the parameters employed in the calculation (which are not 
accurately known) and by the statistical error in lattice simulations. This implies that A (0) 
is determined with a large error (within an order of magnitude for the transition B — > p), an 
uncertainty that heavily affects the analysis of the accuracy of the factorization approach. 
This error is even more important in the study of light SU(3) f breaking effects, an argument 
which has recently prompted a number of interesting investigations [2^ . 

The second point concerns the functional dependence of A on q 2 , which is needed in 
computing nonleptonic transitions. The assumption that, also for low values of q 2 , the 
form factor Aq is dominated by the nearest resonance requires an explicit check, since it is 
known that, for the transitions B,D —>■ p, K*, the form factors Ai and A 2 obtained by QCD 
sum rules appear to be nearly independent of q 2 up to rather large values of the squared 
momentum transferred [15,16|,21 



Both these difficulties can be avoided by a direct calculation of A (q 2 ). This is the 
problem we address in the present paper: we compute A (q 2 ) for the transitions B —>■ p, 
D — > p, K* and B — > D* using three-point function QCD sum rules. In Sect. 2 we derive 
the sum rule for A (q 2 ), and in Sect. 3 we collect our numerical results, together with a 
comparison with the outcome of other calculations. 

The last point we analyze in this note is a proposal, put forward in ref. p3|-f25|, to relate 
the branching ratio of the rare transition B — > K*j to the spectrum of the semileptonic 
decay B — > plv. This connection could have interesting phenomenological consequences, 
e.g. for the measurement of V u b ; it is based on the equality, obtained in the infinite heavy 
quark mass limit, Ti(0) = Aq^ p (0), where T\ is the form factor governing B — > K*j. We 
shall discuss the validity of this relation and show that there is a sizeable deviation to be 
taken into account. 
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II. THE FORM FACTOR A FROM THREE POINT FUNCTION QCD SUM 

RULES. 

In order to calculate the form factor A Q (q 2 ) in eq.([1.8|) using three-point sum rules, let 



us consider the correlator p6|j27 



T^T K \p,p',q) =i 2 \dx dy e^--**) < 0\T{ Ju (x)A,(0)jl(y)}\0 > . (2.1) 

Here the weak current = 57^750 represents the flavour changing axial current in 
(|Q|), while the quark currents j u = d^ v s and j'5 = dij^c interpolate K* and D mesons, 
respectively, and have a non-vanishing matrix element between the vacuum and K* and D 
states: eq.(|1.6|) and 

< 0\j 5 \D(p) >= fD~~~~ , (2.2) 

m c 

where m c is the mass of the charm quark (hereafter we put the masses of the light u and 
d quarks to zero). The product q^T^ u can be decomposed in two independent Lorentz 
structures: 

q» T, u = -i [P v T\(pV 2 , Q 2 ) + q» Up 2 ,p' 2 , q 2 ) } , (2.3) 

with P u = (p+p')v, and the two scalar functions Tj (i — 1, 2) can be represented by a double 
dispersion relation: 

Ti(p 2 , p' 2 , q 2 ) = 1 [ ds ds'-^ 1 — iw 7 ' 9 L + subtractions (2.4) 
yziT) J (s — p z )(s' — p u ) 

with the spectral functions p[ had ^ expressed in terms of physical states: 

p{ had) = pi res) + p{ hr) (25) 

In ( |2.5|) es ^ represent the contribution of the lowest lying resonances (in the zero width 
approximation) and are expressed in terms of the form factor Aq^ k * we are interested in: 

P { r\s,s\q 2 ) = (2n) 2 f K *f D ^(m 2 K * - m 2 D + q 2 )A^ (q 2 )5(s - m 2 D )5(s' - m 2 K .) (2.6) 



pt S \s, s', q 2 ) = (2ir) 2 f K *f D ^°(3m 2 K « + m 2 D - q 2 )A^ K \q 2 )6(s - m 2 D )6(s' - m K 



2 

D 

(2.7) 



while p\ r ^ take contribution from higher states and from the hadronic continuum. 



The correlator (2A) can also be calculated in QCD by the operator product expansion 
(OPE) in the region of large Euclidean momenta p 2 , p' 2 . In this expansion the most singular 
term is represented by the perturbative term, which can be obtained evaluating a triangle 
diagram. It has a dispersive representation similar to 



5 



T (pert) (v 2 v' 2 a 2 ) = 1 { ds ds' P^^^'^ 2 } . (2 8) 

the explicit form of the spectral functions p\ pert ^ can be found in the appendix. The next 
terms in the OPE can be read as: 

T M = J2d ( r\p 2 ,p' 2 , q 2 )<0\O n \0> (2.9) 

n 

where the local operators O n , written in terms of quark and gluon fields, are ordered accord- 
ing to their increasing dimension; their vacuum matrix elements (condensates) parametrize 
the effects of the nonperturbative QCD vacuum. Here we only consider the contribution of 
the operators of dimension D = 3 and D = 5: 

T/ np) = 4 3) <qq> +d? ] < qg s oGq > ; (2.10) 

< qq > and < qg s aGq > are the quark and the mixed quark-gluon condensates, respectively; 
the coefficients df^ and df^ can be calculated in perturbation theory. 

The hadronic and the QCD representations of the correlator (|2.1| ) can be used to derive 
two sum rules for A (q 2 ). As a matter of fact, we invoke quark-hadron duality and assume 
the equality: 

{hr) (s s' a 2 ) r o {pert) (s s' a 2 ) 

D' (s — p 2 ){s' — p' 2 ) JD' (s — p 2 )(s' — p' 2 ) 

in a region D' above some threshold Sq, s' . This assumption can be realized if we adopt the 
model: 

p? r \s,s',q 2 ) = pf eH \s,s',q 2 )(l - Q(s - s)Q(s' - *')) • (2.12) 
Using ( 2.12|) the following sum rules for A (q 2 ) can be written: 

o {res) (s s' a 2 ) 1 r o {pert) (s s' a 2 ^ 

ds ds' Pi } : ,q L = / ds ds' A }: ,q 



(27t) 2 Jd (s - p 2 )(s' - p 12 ) (27t) 2 Jd (s - p 2 )(s' - p' 2 ) 

+ df ] <qq> +df ] < qg s aGq > (2.13) 

where the region D is fixed by: 

.2 



m 2 < s < s (2.14) 



and 



m 2 s < s'< s' 



si(s,q 2 ) < s'< s' + (s,q 2 ) (2.15) 
with (m s is the strange quark mass): 
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s — m. 





1/2' 



(2.16) 



In ( |2.13| ) we have omitted the subtraction terms, that can still be present. They are removed 
by performing a double Borel transform to both sides of ( 2.13 ) in the variables P 2 = —p 2 



and P' 2 = -p' 2 : 



n+l 



_W_p2)n+l 

n\ K ' \ dP 2 



1 



d 



m+1 



mV ' \ dP' 2 



(2.17) 



in the limit P 2 ,P /2 — > oo, n, m — > oo, M 2 = P 2 /n M' 2 = P' 2 /m fixed. The resulting sum 



rules for Aq^ k * read: 



H 



(2vr) 2 Jd 



ds ds' pf ert ' (s, s', q 2 ) exp( 



Is — m 



D) 



is' — m 2 



K* 



M 2 



+ 



A" 



rf =A <qq> +rr =b < qgs^G^—q > 



cxp 



M' 2 



M 2 



.m* D — m 2 m 2 K * — m" a . 



M' 2 



(2.18) 



(i = 1,2); the coefficients Hi, Tf =3 and rf =5 are collected in the appendix. 

Eq. (|2.18f) represents two independent sum rules for Aq^ k * (q 2 ), that will be analyzed in 
the next section. Now, before discussing the parameters and the criteria we have used to 



compute Aq~* k * , let us consider the t dependence in ( |2.18| ). In principle, the operator prod- 



JJ) in the region 
we can also consider positive 



uct expansion can be reliably applied to the evaluation of the correlator 
af large euclidean values of q 2 
values of q 2 provided that the occurrence of non-Landau singularities either is avoided or 
is carefully taken into account. Such singularities remain far from the integration region of 
the dispersive integral if q 2 is small with respect to the physical threshold in the t channel: 
Uh = ( m c + m s 
and positive q 2 

Taking into account the above considerations, we study the form factor A using eq. fl2.18Q 
for small positive values of q 2 , and not in the whole physically accessible t-region. However, 
the information we obtain on the t— dependence of the form factor is enough to describe a 
large number of nonleptonic D (B)— meson decays. 



2 ; their presence is shown up by the appearance of branch points (for large 
in the spectral functions p 4 - per ^. 
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III. NUMERICAL ANALYSIS OF THE SUM RULES 



The numerical analysis of the sum rules ( |2.18| ) is performed using the following values 



for the quark condensates, taken at a low renormalization scale (yU — 1 GeV) |2 

< q q > = (-230 MeVf 
\ a 

< qgs^G^—q > =m 2 < qq > (3.1) 

with wig = 0.8 GeV 2 . The condensates can be evaluated at higher scales using the leading- 
log approximation for the anomalous dimension of the quark and of the mixed quark-gluon 
condensates; however, this rescaling does not affect sensitively the numerical result for A . 

As for the quark masses, we use: m c = 1.35 GeV and m s = 0.16 GeV; moreover, we 
use rrib = 4.6 GeV in the calculation of form factor connected with the transition B — > p. 
The leptonic constant fx* can be derived from the measurement of the branching ratio 
t -> K*u r : f K * = 0.22 ± 0.01 GeV; we use f D = 195 ± 20 MeV for the leptonic constant 
of the D meson. 

A comment on this value of fu is in order. This value comes from two-point QCD sum 



rules, including radiative corrections at the order 0(a s ) P5 , p9| ; these corrections are at level 
of 13 — 15 % for f D , f D * and play a major role in f B ; including the corrections the results are 
in agreement with lattice QCD. On the other hand, the QCD expression of the three-point 
correlator has been computed at zero order in a s . Several authors 15H , when computing 



the form factors, adopt the strategy of taking the ratio of three and two-point functions, 
calculated at the same order in a s in order to be consistent, trying to reduce the uncertainty 
deriving from higher order contributions. However, the sign of the radiative corrections to 
three point functions is not known on general grounds; for this reason our attitude is to use 
the best known values of the leptonic constants, e.g. the values on which QCD sum rules 
and lattice QCD agree. 

The effective threshold Sq and s' must be chosen in a range of values between the mass 
squared of the lowest lying resonances and the first excited states. They have been fixed 
by studying two-point sum rules, for the calculation of static properties as the mass of the 
leptonic constant of the particles, and three-point functions sum rules in the calculation of 
the semileptonic D — > K*lv decay [15]: sq = 7 — 8 GeV 2 for the channel of D meson, and 



s' = 1.5 — 1.7 GeV 2 for the channel of K*. The variation of the values of sq, s' induces a 
variation of the predicted form factor A that can be considered a theoretical uncertainty in 
the final result. 

Using the parameters above, A®^ K * (q 2 ) can be obtained from ( |2.18| ) as a function of 
the Borel parameters M 2 and M' 2 . However, since these variables are unphysical, we search 
a region in M 2 , M' 2 where Aq does not depend on them (stability plateau). In this region 
other conditions must be verified. The first constraint consists in checking a hierarchical 
structure of the various terms of ( |2.18| ): 

rp(pert) > j,(D=3) > rp(D=5) . 2) 

only if the OPE displays this structure we can hope that higher order power corrections in 
( |2.18| ) can be safely neglected. 



S 



Another condition is connected with the approximation ( |2.11|) . Since we are not guar- 
anteed that quark-hadron duality starts already in correspondence to the first hadronic 
excitations, we must choose M 2 , M' 2 small enough to enhance the contribution of the lowest 
lying states in the sum rule, and to exponentially suppress the contribution of the continuum. 
We check this condition by verifying that 



D 



ds ds'p 



(pert) 



s,s',q 2 ) 





(s-m 2 D ) 


(s'-™ 2 K *) 


> a 




e 


M 2 


M< 2 


L 













ds ds'p 



(pert) , 



4)_ 
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(3.3) 



where D\ is an integration region larger than D (it extends to 10 so and 10 s' ) and a = 
0.4 — 0.5. From the condition ( |3.3| ) an upper bound to M 2 , M' 2 can be determined. 

In principle both the sum rules (i = 1, 2) in Q2.18D can be used to calculate A . However, 
it turns out that only for the first sum rule (i = 1) the condition (|3.3|) is fulfilled. In the 



second sum rule the constraint (|3.3| ) is never verified as a consequence, probably, of the 
major role played by higher states. For this reason we exclude the second sum rule in our 
analysis and use only the first one to determine Aq. 

Let us consider the form factor A®^ K (q 2 ). Our results, obtained for M 2 = 3 GeV 2 and 
M' 2 = 1.5 GeV 2 , are depicted in fig.(l) for several sets of parameters. At q 2 = we have: 
Aq~* k *(Q) = 0.58 ±0.05, where the uncertainty comes from the variation of the parameters. 
For positive values of q 2 we observe an increasing of the form factor; if we use the simple 
pole formula: 



aD- 
^0 



^0 



>K* 



(0) 



1 — q 2 /m P 



(3.4) 



to fit the q 2 dependence, we obtain that the fitted pole mass is: trip = 1.65 GeV, which is 
not far from the experimental mass of the particle having the same quantum numbers of the 
pole: m Ds = 1.969 GeV. 

Let us now consider the transition D —>■ p. The constant /„+, computed from the decay 



p u — > e + e and assuming the isospin symmetry, is given by: f p + =216 MeV. Using 
s' Q = 1.3-1.5 GeV 2 we get: A^ p (0) = 0.52 ±0.05. The fitted pole mass is m P = 1.6 GeV, 
to be compared with wid = 1.865 GeV which is the mass of the nearest pole in the t channel 
(fig.(2)). 

From the above results we would get: r = A®^ K * (0)/A ^ p (0) = 1.12 ± 0.11; however, 
the light SU(3) p breaking effects can be better estimated by studying the ratio of the sum 
rules for the transitions D — > K* and D — ► p, which is stable with respect to the variation 
of the input parameters. We get r = 1.10 ± 0.05, i.e. an SU(3) P breaking effect at the level 
of 10%. 

We have also computed the form factor Aq^ p appearing in nonleptonic B meson decays 
which are interesting due to their connection to the measurement of V u b, using so = 33 — 
36 GeV 2 and f B = 180 ± 20 MeV as obtained by two point function QCD sum rules, we get 
Aq^ p (0) = 0.24 ± 0.02 and a pole mass mp = 5 GeV, to be compared to mg = 5.275 GeV; 
the form factor is depicted in fig. (3) in the range q 2 = — 10 GeV 2 . 

Finally, for the heavy-to-heavy meson transition B — » D* we get (using fp,* = 250 MeV 
Pi): A o^ D *(°) = 0.65 ±0.05. 
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Let us now compare our results with existing predictions for Aq(0), obtained from the 
values of Ai(0) and A 2 (0) in the corresponding channels (see Table I). The first observation is 
that our predictions appear to be smaller than from potential models |7|, [il||T4|| (the exception 
is Aq~* p in [13[]). On the other hand, the SU(3) breaking effects are of the same size as in 
our approach. 

The comparison with QCD sum rules predictions |]lT>l , |T6l shows the problem we have 
mentioned in the introduction: the uncertainties in the results obtained using ^i(O) and 
A 2 (0) are so large that they obscure the real value of A (0), whereas the error in our 
predictions is at the level of 10%, and therefore a substantial improvement has been obtained. 

It is interesting to compare our predictions to the outcome of [|19| obtained in a framework 
based on heavy-quark and chiral symmetries, using as an input the experimental data on 
D — > K*£u. Although there is a remarkable agreement on D — > K* and B — > p, the light 
SU(3)p breaking corrections connecting D — > K* to D — > p are of similar size than in our 
approach but with opposite sign. 

The last point we would like to mention is that, in our calculation, the form factors A in 
the different channels increase for positive values of q 2 , and that their functional dependence 
is compatible with the simple pole behaviour dominated by the nearest resonance in the t 
channel. 



IV. B K*7 VERSUS B -> piv 

In this section we discuss a procedure, proposed in refs. |23| , |24|| , to relate the width of the 
radiative B — > K*j transition to the spectrum of the Cabibbo suppressed semileptonic decay 
B — > piv. This relation can be used to constrain the electroweak parameters involved in 
these processes and to reduce their dependence on the models for the hadronic form factors. 

The radiative rare B meson decays, like B — > K*j (K* = K*(890), Ki(1400), etc.), have 
been extensively studied, from the theoretical standpoint, since they have a peculiar role in 
the precision tests of the quark sector in the Standard Model (SM) and probe the effects of 
new physics beyond the Standard Model |j30| . Within SM they are induced by the one-loop 
electromagnetic penguin operator b — > 57 (for m s « m b ) f3T| : 

Heff(b -> s) = ^Cm b e»sa^J-^q v b (4.1) 

where the constant C contains the dependence on the Cabibbo-Kobayashi-Maskawa matrix 
elements, on the QCD correcting terms and on the ratio m top /mw (the explicit formulae 
can be found in ref. . 

The description of the exclusive decays requires the knowledge of form factors. For 
example, the relevant matrix element for the transition B — > K* (890)7 can be written as 
follows: 
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(K*(p',\)\sa pu q 



,(1 + 75) 



b\B{p)) =ie„ vprj e* v p»p"> T^q 2 ) 



+ 



m, 



(p 



(m 



B 



m 2 K *) 



Ap+p') 



Uq 2 ) . (4.2) 



The form factors Tj(g 2 ) require a nonperturbative evaluation. However, in the static limit for 
the b quark and using the SU(3)p symmetry for the light quarks, a number of interesting 
relations can be derived between Tj and the form factors V and A{ governing B — > plv 
23| , |32[| . Although such relations are obtained at zero recoil {q^ax 



[m B - m K *, P 

-2 



it has 



been argued that they continue to hold in the whole physically accessible q range p3| , |33| , |3l 
and therefore that one is able to relate the width of B 
B — > plv. An example is represented by the ratio: 



if* 7 to observable quantities in 



Y{B -> K*i) 



lim 



1 dT{B -> pev) 



-1 



. <7 — >0, curve q 



dE p dE e 



4tt 2 


H 


2 {m\ — m 2 K *) 3 


^2 1 


v ub 


2 m B 



(4.3) 



that Burdman and Donoghue [23|] predict to be independent of any form factor. In ( |4.3| ) 
curve means the region in the Dalitz plot where q 2 = AE e (mB — E p — E e ); the factor 
r) includes the QCD corrections to the b — > S7 decay and the other relevant electroweak 
parameters jp3"| . 

Eq. ( |4.3| ) is interesting, both from the experimental and the theoretical point of view, since 
it would provide, for example, a method to measure V u b without referring to calculations of 
the form factors, from the measurement BR(B — > K*j) = (4.5 ± 1.0 ± 0.9) x 10~ 5 [^S] and 
the spectrum of B —>■ plv. The drawback, pointed out by O'Donnell and Tung pi |, is that 
the semileptonic decay spectrum vanishes at q 2 = on this curve, with the consequence that 
the number of experimental points required to measure the denominator in ( |4.3|) could be 
insufficient for a reliable determination of V u b. 

A more convenient observable [^,0 is represented by the ratio: 



R(B -> K*j) 



dT(B -> pW t ) 



dq 2 



1927T 3 



[mi 



m 



K" 



ml 



r<2 



Wbul 2 



(m 2 B — m^) 3 



(m 2 — m 2 s ) z 



ixr 



(4.4) 



In (|4.4j) R(B — > K*^) = T(B — > K*^)/T(b — > 57) is independent of the electroweak param- 
eters appearing in the effective Hamiltonian, and only depends on Ti(0): 



R(B -> K*~f) 



ml(m 2 B — m 2 ^*) 3 
m 3 B (ml — m 2 s ) z 



T!(0)| : 



having used the relation T 2 (0) = |Ti(0). The factor X is given by p5| 



1 



(m B + m p ) Ti(0) 



(m B + m K *) Aq^ p (0) 



(4.5) 



(4.6) 



The main reason for studying the ratio ( [4.4|) is that dT(B — > plvi)jdq 2 in ( [4.4| ) does not 
vanish at q 2 = 0. Moreover, the factor I should be close to one. As a matter of fact, in the 
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framework of the Bauer-Stech-Wirbel model M the value X = 1.12 can be derived, whereas 



the constituent quark model [37] provides the value X = 1.09 — 1.18. The relation X = 1 
has been obtained by the authors in ref. |24| by considering the transition bq — > Qq in the 
limit where both the b and Q quarks are heavy, applying spin symmetry relations; then, it 
has been argued that this relation still holds for a light Q quark in the weak binding limit 
for the meson. 

The factor X can be derived by three point QCD sum rules, using realistic values of the 
quark masses in a fully relativistic approach. Using the determination T 1 1 (0) = 0.35 ± 0.05 
|3"3| and the value Aq^ p (0) derived in the previous section, we get: X = 1.43 ±0.32 . Notice 
that, in evaluating the error for X, we have taken into account the correlated errors in the 
numerator and in the denominator of X. 

The uncertainty can be reduced by studying the ratio of the sum rules determining 
7i(0) (the formulae can be found in ||33|| ) and Aq^ p (0). We find a wide stability window in 
correspondence to the result: X = 1.3 ± 0.1 . 

Therefore, the QCD sum rules prediction differs by 10 — 30% from the values derived 
in [p^ , pS| ; this deviation still is model dependent, in the sense that it comes out from an 
explicit QCD sum rules calculation and not from symmetry arguments, and therefore it 
should be checked in a different approach, e.g. by lattice QCD. 

In any case, we feel that, modulo this uncertainty on the value of X, the idea of relating 
B — > K*~f to B — > ptv provides us with an interesting alternative method to measure V u b. 



V. CONCLUSIONS 

We have applied three-point function QCD sum rules to calculate the form factor A 
appearing in the matrix element of the flavour changing axial current between the D and B 
mesons and a vector meson state. This form factor plays an important role in the calculation 
of non leptonic B and D decays to PV states in the hypothesis of factorization. 

The explicit calculation has allowed us to determine A Q at zero momentum transferred, 
with better accuracy than in the indirect determination from A\ and A 2 . We have also com- 
puted the q 2 dependence, finding an increase for positive values of the squared momentum 
transferred, compatible with a simple polar behaviour. 

Finally, we have considered the role of Aq in relating the rare B — > K*^ decay with the 
semileptonic B — > piv spectrum, in the search of methods for measuring the matrix element 
V ub . 
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APPENDIX A: PERTURB ATIVE AND NON-PERTURBATIVE TERMS. 



The perturbative spectral functions pf (s, s', q 2 ) in fl2.8|) , at the lowest order in a s , can 
be obtained by the Cutkosky rules applied to a triangle diagram: 



PiO, s r , q 2 ) = -^={(m c + m s ) (A' - A) 



2(A s' + A' s) —u (A' + A) 
A 



x 



(to c + m s )(s — s' + q 2 ) + 2m c {m 2 c — to^) | 



P2 



(s, s', q 2 ) = ^{K + m s ) (A' + A) + 



2(A s' — A' s) + w (A — A') 
A 



x 



(m c + m s )(—s + s' — q 2 ) + 2m c (m 2 — to 2 ) } 



where A = s — to 2 , A' = s' — m 2 , u = s + s' — q 2 and A = m 2 — 4 s s'. 

The Borel transformed coefficients of the non-perturbative D = 3 terms read: 



n 



(D= 3 ) _ mg - m 2 



(Al) 



,(D=3) _ (™ C + ™ S )' 



the coefficients of the condensate of dimension D = 5 terms are: 

r (D= 5 ) = mg(m 2 -mg) m 2 (m 2 - m 2 c ) (m 2 c - m 2 ) (mg - to]) 
1 8M 4 8M' 4 6M 2 12 M' 2 

(m 2 — mg) (2m 2 — mg) — q 2 {m s + m c )(2m s — m c ) 



+ 



12 M 2 M' 2 



and 



r (£>=5) 
1 2 



(A2) 



mg(m,, + m c ) 2 mg(m s + m c ) 2 (mg + 3m 2 + 4m s m c ) (3mg + mg + 4m s m c ) 



+ 



8M 4 + 8M' 4 12 M' 2 

(m 2 + mg — q 2 ) (2m 2 + to 2 ) + m s m c (Am 2 + 2m 2 — 3g 2 ) 



6M 2 



12 M 2 M' 2 

Finally, the coefficients Hi in ( 2.18|) read: 



Hi 



2m c 



and 



Ho 



fK*fDm 2 D (m 2 K , -m 2 D + q 2 ) 



2m r 



(A3) 



f K *fDm 2 D (3m 2 K * +m 2 D - q 2 ) ' 
The formulae for the transitions D — > p, i? — > p and -B — > D* can be easily obtained. 



(A4) 
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TABLES 



TABLE I. The form factor Aq(0) of various transitions in different models. 



Transition 







@0I 



This paper 



D ^ K* 

D^p 

B^p 



0.74 
0.68 
0.28 



0.91 



0.8 

0.85 

0.14 



0.45 ± 0.30 
0.57 ±0.40 
0.79 ± 0.80 



0.77 ±0.29 



-0.57 ±0.65 



0.59 
0.74 
0.24 



0.58 ± 0.05 
0.52 ±0.05 
0.24 ±0.02 
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FIGURE CAPTIONS 



Fig. 1 

The form factor A^ K *(q 2 ). The curves refer to different sets of parameters: s = 7 GeV 2 
and s' Q = 1.5 GeV 2 (continuous line), sq = 7 GeV 2 and s' = 1.7 GeV 2 (dashed line), 
So = 8 GeV 2 and s' = 1.5 GeV 2 (dotted line), s = 8 GeV 2 and s = 1.7 GeV 2 (dashed- 
dotted line). M 2 = 3 GeV" 2 , M' 2 = 1.5 GeV 2 . 



Fig " 2 

The form factor A ^ p (g 2 ). The curves refer to the sets of parameters: s = 7 GeV 2 and 
s{, = 1.3 GeV 2 (continuous line), s = 7 GeV 2 and s' = 1.5 GeV 2 (dashed line), s = 8 GeV 2 
and s' = 1.3 GeV 2 (dotted line), s = 8 GeV 2 and s = 1.5 GeV 2 (dashed-dotted line). 
M 2 = 3 GeV 2 , M' 2 = 1.5 GeV 2 . 



Fig. 3 

The form factor A ~* p (g 2 ). The curves refer to the sets of parameters: s = 33 GeV 2 
and s' = 1.3 GeV 2 (continuous line), s = 33 GeV 2 and s = 1.5 GeV 2 (dashed line), 
s = 36 GeV 2 and s' = 1.3 GeV 2 (dotted line), s = 36 GeV 2 and s' = 1.5 GeV 2 (dashed- 
dotted line). The values of the Borel parameters are: M 2 = 8 GeV 2 , M' 2 = 2 GeV 2 . 
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